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We predict the polarization of cosmic microwave background (CMB) photons that results from a 
cosmic bubble collision. The polarization is purely -E-mode, symmetric around the axis pointing 
towards the collision bubble, and has several salient features in its radial dependence that can help 
distinguish it from a more conventional explanation for unusually cold or hot features in the CMB 
sky. The anomalous "cold spot" detected by the Wilkinson Microwave Anisotropy Probe (WMAP) 
satellite is a candidate for a feature produced by such a collision, and the Planck satellite and other 
proposed surveys will measure the polarization on it in the near future. The detection of such a 
collision would provide compelling evidence for the string theory landscape. 
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Recently, considerable attention has been focused on the question of the observability of cosmic bubble collisions. 
These collisions occur in models in which our entire Hubble patch — in fact, the entire region of the Universe in our 
vicinity describable by an approximately homogeneous and isotropic metric — exists inside a bubble created by a first- 
order phase transition from an eternally inflating false vacuum, and our bubble is struck by another that nucleated 
(from our parent false vacuum) nearby. They are of interest in part because their existence is predicted by the string 
theory landscape [H |2] and is intimately related to its solution of the cosmological constant (CC) problem. 

In the landscape the CC problem is solved by anthropic selection; that is, we find ourselves in a region of small 
vacuum energy because such phases exist in some places in the landscape, and it is only inside bubbles of these 
phases that structures like stars and galaxies can form ^3^. One expects such small CC bubbles to be created by 
nucleation from a rapidly inflating false vacuum. If so, and if the instanton mediating the decay has the symmetries 
of the Coleman-de Luccia solution [4|, the inside of the bubble is homogeneous and isotropic with negative spatial 
curvature (it contains inside it a complete, spatially infinite open Friedman-Robertson- Walker Universe). Viewed 
from the outside the bubble expands at nearly the speed of light, but remains immersed in a bath of even more 
rapidly inflating false vacuum. This phase is necessarily metastable and will occasionally decay to form other bubbles. 
When another bubble — either of the same type or another — nucleates within one false-vacuum Hubble length of the 
wall of our bubble it collides with it, leaving an interesting and highly characteristic imprint on various large-scale 
cosmological observables [H [3] . 

Modulo the (substantial) theoretical uncertainties inherent in the above paradigm, the existence of such collisions 
is a prediction of string theory. Hence the possibility of observing such an event presents an opportunity to greatly 
enhance our understanding of both cosmology on the largest possible length scales and of fundamental microphysics — 
and as such, it should be pursued vigorously. The results are of interest outside of the context of string cosmology 
as well because the existence of such collisions is generic in any theory with slow-first order phase transitions coupled 
to gravity. In addition, the techniques developed to analyze these collision spacetimes are of utility for the study of 
anisotropic and inhomogeneous cosmologies more generally. 



In general, the expected number of potentially observable bubble collisions with bubbles of type j was estimated 
in [7]; it is 
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Symbols and Definitions 


To 


Average photon temperature today (~ 2.726 K) 


Tdc{x,y,w) 


Photon temperature at decouphng (t — tdc) 




Comoving distance from Earth to the decouphng/recombination surface 




Comoving distance to a scattering electron 


a 


The ratio D^/ Ddc 


Xc 


Comoving location of the collision lightcone at decoupling 


X 


Linear slope of the temperature perturbation at decoupling 


A 


Slope of the "effective" perturbation including Sachs- Wolfe (see text) 


Oc 


Angular radius of the affected disk on the CMB temperature map 




cos 6 



where 7^- is the decay rate of the false vacuum to vacua of type j (in units of the false vacuum Hubble rate) , Vf is 
the parent false vacuum energy density, Vi is the energy density during slow-roll inflation in our bubble, and VL^ is 
the magnitude of the negative spatial curvature in our bubble today. This is the number of collisions with future 
lightcones that bisect the part of the last scattering surface visible in the Cosmic Microwave Background (CMB) sky 
today (collisions whose future lightcones have not yet affected the CMB are unobservable, and those that encompass 
our entire sky and more are probably hard to detect). This number could be significantly greater than one without 
fine-tuning (although it appears unlikely to be of order one). 

Among the theoretical uncertainties is our ignorance of the decay rates of our parent false vacuum. If all the rates 
are extremely slow, the spacetime region in our bubble that we can observe today (despite its age of 13.7 Gyr) may 
not yet have been affected by a collision with another bubble. Another caveat is the amount of slow-roll inflation in 
our bubble. After N e-folds of slow-roll inflation, the spatial curvature is rjfc ~ e^^^ where N* is 60 (up to 
logarithmic dependence on factors like the reheating temperature (see e.g. [S]), so increasing the number of e-folds 
of inflation by AA^ decreases \/Vlk by a factor of e~'^^. Inflation also reduces the magnitude of the effects even 
for bubbles that bisect the CMB sky (by decreasing the brightness of the perturbed disk they create, for example). 
Hence large amounts of inflation both increase the distance between collisions by expanding our bubble and reduce 
the magnitude of their effects by inflating away their remnants. Therefore, like any pre-inflationary relic, bubble 
collisions become exponentially harder to observe as the number of e-folds of inflation increases. 

However, it appears that large numbers of e-folds of inflation require fine-tuning to attain in string theory, and 
at the same time that there is an anthropic requirement that there be sufficient inflationary expansion to solve the 
flatness problem [5]. One then expects that the number of e-folds of inflation was not too much larger than the 
required minimum. The theoretical uncertainties are too great to say anything more deflnite, but it does not appear 
that a large degree of flne-tuning is required for the effects of these collisions to be observable. 

B. Previous Work 

The cosmological implications of bubble formation and collision were studied recently in [5H71 [SHTT] . Speciflcally, 
[5] derived the approximate CMB temperature anisotropy due to a collision — the effect is conflned to a disk that is 
either hot or cold relative to the average, with an intensity that decreases linearly with cosine of the angular radius 
from the center of the disk, reaching zero at its edge. That result will be our starting point here, but rather than 
reproducing the analysis of [5] we will begin with the collision temperature perturbation at decoupling and study its 
effects on CMB polarization. The effect of the collision on the peculiar velocities of large scale structures was analyzed 
in [6j. Previous work on CMB polarization and the WMAP cold spot appeared in [T5] . 

C. Summary of Results 

We compute the contribution to CMB polarization from a cosmic bubble collision, treating the collision's effects 
using a technique based on [5 . We find that the polarization is purely -E-mode (as expected for a scalar perturbation) , 
and is radial (azimuthal) for collisions that produce cold (hot) disks in the temperature map. If the polarization 
resulted entirely from Thomson scattering off electrons at a single value of the redshift z, the effects would be 
confined to an annulus containing the edge of the affected disk in the temperature map. With a realistic model for 
scattering, the polarization is largest near the edge of the temperature disk, and drops to zero at its center and as 
one moves radially out from the edge. 
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Our result is accurate for collision disks larger than a few degrees, but cannot be relied on quantitatively for features 
at sub-degree scales (see Sec. |v]for details). If the WMAP cold spot [I9H23] is the result of a bubble collision, the 
corresponding polarization pattern is strong enough to be detected by the Planck satellite. We do not analyze the 
question of with what confidence one could rule out a random Gaussian origin were the measured polarization to be 
consistent with our prediction. However, we do emphasize that the temperature and polarization pattern we predict 
has a striking planar symmetry in three dimensions (inherited from the physics of the collision) which may be difhcult 
to mimic with random fields. 



II. BUBBLE COLLISION BASICS 



We focus on a collision scenario where the bubble we inhabit formed by a Coleman-de Luccia transition from an 
eternally inflating parent vacuum, was struck by a bubble of a third type, and underwent a period of slow-roll inflation 
soon after its formation. This is schematically shown in Fig. [T] We assume that the pressures on the domain wall 
formed between the two colliding bubbles are such that the wall accelerates away from our bubble. This assumption 
is valid for collisions between two de Sitter (dS) bubbles if ours has a smaller vacuum energy, and in collisions with 
anti-de Sitter bubbles when certain conditions involving the tension of the wall and the vacuum energies are satisfied 
[TTJ [ini mi m] . WhUe more general scenarios are possible, these assumptions can lead to a cosmology consistent with 
current observations. Finally, we assume that the collision lightcone divides the part of the decoupling surface that 
we can see in the CMB sky. 

We wish to investigate the polarization of CMB photons induced by the collision. A complete treatment would 
include the non-linear evolution of the cosmology from the highly anisotropic and inhomogeneous initial conditions 
set up by the collision. However, since observations constrain the magnitude of the collision signal to be small in the 
region of the universe we can observe today, perturbation theory is sufficient to characterize its effects on most (if not 
all) signals. The precise formulation of cosmological perturbation theory in this context and a detailed analysis of the 
evolution of cosmological perturbations in such a bubble collision Universe will be discussed in future work 

In this paper we instead make several approximations so that an analytic analysis that captures the key physics is 
possible. We will treat the effects of the collision starting from the time t = tdc of decoupling between photons and 
atoms (also known as the time of recombination). We assume that inflation lasted long enough so that the negative 
spatial curvature can be ignored, so that to leading order the background metric is approximately 

ds^ = -dt^ +a^{t)dx^. (2) 

Here x = {x, y, w) is a comoving coordinate, the x-axis points towards the collision, and the lightcone of the collision 
event is a plane extended in the y and w directions, and null in the x and t directions (Fig. [T]). Without loss of 
generality the coordinate position of the Earth today is taken to be x = (0, 0, 0) and t — to. 

To model the effects of the collision, we take the perturbation to the temperature of the Universe due to the collision 
to be 

T<j,(x) = (1 + Zdc)To [l + X{x - xc)e{x - xc)] (3) 

at t — tdc- Here Q{x) is a step function, the comoving coordinate of the collision lightcone at decoupling 

(see Fig. [l]), and A is a parameter that controls the intensity of the affected disk in the CMB sky. The temperature 
perturbation is assumed to be adiabatic with no isocurvature contribution. The temperature in the part of the 
Universe unaffected by the collision is T^c = 7o(l + Zdc), where Tq ~ 2.726 K and Zdc — 1089. As we demonstrate 
below, this temperature profile reproduces the results of [S] for the CMB temperature perturbation due to a bubble 
collision. 

So long as the perturbations are small, they may be treated linearly, and (as in 15 J we will assume that the random 
fluctuations produced during slow-roll inflation remain Gaussian and uncorrelated with the collision perturbation at 
lowest order. As such, we will not discuss the usual random contribution to anisotropics but they will also be present 
with their usual amplitude. 

Note that for large x — Xc the effects of the collision become large and cannot be treated perturbatively. This is 
due to the fact that near the domain wall separating the bubbles, the spacetime is strongly affected and cannot be 
approximated by a perturbed FRW metric. However, consistency with the observed isotropy of the CMB requires 
that the perturbation be small in the part of the decoupling surface that falls within our past lightcone today, and as 
long as that remains true our approximation is valid. 

The ansatz adopted here differs somewhat from that of O [B]. Roughly speaking, [S] |2] treated reheating and 
decoupling as surfaces of constant temperature but non-constant FRW time, while here we treat them as surfaces of 
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Bubble Walls 



FIG. 1. Sketch of the causal structure of the spacetime in the bubble collision scenario. The region to the right of the radiation 
shock is within the forward lightcone of the collision and affected by it. The figure shows the trajectories of two photons that 
last scattered at decoupling, and one that last scattered at reionization. The figure is not to scale. 

constant time but non-constant temperature (and gravitational potential). To linear order in the perturbation these 
approaches are very similar (if not simply gauge-equivalent). The end result of all these analyses is a perturbation 
with a spatial dependence that is dictated by the symmetries of the bubble collision and two parameters (the distance 
from the Earth's comoving location a; = to the collision lightcone x — x^., and the slope of the effect inside the 
lightcone A). The precise relationship between these approximations and their regimes of validity will be clarified in 
detail in [5S]. 

A. Temperature Anisotropy 

As a first step, we compute the CMB temperature anisotropy due to the collision perturbation. The temperature 
anisotropy ([S]) gives rise to a gravitational potential perturbation $, and the CMB temperature observed in a given 
direction today is the sum of the intrinsic temperature and gravitational potential at the corresponding point on 
the sphere defined by the earth's past lightcone at t = tdc- For an adiabatic perturbation the gravitational potential 
•^dc - i~2/3){ST/T)dc [261E2 , so the temperature anisotropy observed today is 6T/T ~ {-1/2){6T/T)dc^ Therefore 
for our purposes we can incorporate the Sachs- Wolfe effect by defining an "effective temperature" anisotropy based 
on ([3f: 

T{9,^) = Tdc^eff{Ddcn)/{l + Zdc)=To[l + XDdc{fi-fic)eU^~^ic)], (4) 

where A ~ —A/2 and 

rrfc,e//(x)/(l + Zdc) = To [1 + ~ X,)eix - X,)] (5) 

is the effective temperature perturbation taking the gravitational potential perturbation into account, Ddc = 
/;2 dt'/ait') is the comoving radius where the past lightcone of the Earth today intersects the decoupling sur- 
face, n is a unit vector pointing from the Earth in the direction (0,0), /i — cos 6, and /ic = cos 0c = Xc/Ddc is the 
angular radius of the collision disk in the temperature map. We have chosen coordinates so that the collision disk is 
centered at = (/i = 1). 

From Q we see that the effect of the collision is to create a disk on the CMB temperature map, inside of which 
the temperature depends linearly on the cosine of the angle from the center, and outside of which there is no effect. 



^ In general there are additional contributions from velocity perturbations, but for features larger than about a degree these are subdom- 
inant and we will ignore them. 
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The angular radius of that disk is 6c = cos~^ {Xc/Ddc), and the maximum temperature deviation is (AT/T)max = 
^{Ddc — Xc)- The disk can either be hot or cold depending on the sign of the parameter A; in terms of the microphysics 
of the collision this depends on the characteristics of the other bubble, the domain wall, and the inflaton (see [S] for 
a discussion). With the mapping jic — >■ —xt, ^{Ddc — Xc) — > M — 1 (the collision disk in this paper is centered at 
9 = 0) the temperature anisotropy Q coincides with the one discussed in [5^ . 



III. POLARIZATION BASICS 



The CMB is polarized because of the Thomson scattering of CMB photons off electrons [28l [29] . This scattering 
occurs primarily at redshifts centered at decoupling (z^c ^ 1100) and reionization {Zre ~ 10) as shown in Fig. [T] 
The contribution from decoupling depends primarily on known atomic physics |30j . but there is currently significant 
uncertainty in the reionization history. We will derive an expression for the contribution to the polarization from the 
collision in terms of a general visibility function, and then evaluate it using two fiducial reionization models. 

As discussed above in Sec. |II A| a bubble collision introduces an inhomogeneity in the temperature at the decoupling 
surface. As a result, some of the electrons in our past lightcone see a quadrupole temperature anisotropy from their 
CMB skies. Photons scattered off such electrons then give rise to polarization of the CMB. This section briefly reviews 
the physics of this effect ^5112^ 1 ^ 15^. 

Consider a typical electron in the Universe at some time t = and comoving location Xg. Photons scattered by this 
electron and observed by us today will be polarized if they originate from a temperature distribution such that the 
anisotropy pattern in the CMB sky seen by the electron has a non-zero quadrupole moment. Therefore as a first step, 
we should compute the quadrupole moment of the temperature anisotropy on the intersection of the past lightcone 
of the electron with the decoupling surface. For an electron at comoving distance De, that intersection is a sphere of 
comoving radius Re = Ddc — (see Fig. [2]). Therefore for an electron at position Xei 

2,„(x,) = j dn, r;,„(n,) (ST/T)^ = ^ ^^^^^^ J dh^ Y*^(ne) Tdc.effi^. + i?efie) (6) 

Here fig is a unit vector pointing from the electron at position Xg in the direction (0e, 0e), and Tdc,e//(xe + Re^e) is 
the effective temperature distribution incident on the electron (c./. ([5|). We have used the independence of 6T/T on 
redshift and the orthogonality of the monopole and quadrupole in the second equality. 

Photons scattered off the electron carry linear polarization proportional to the quadrupole of ([6|. To find the 
polarization at a fixed direction n on our sky, one must integrate the contributions of all electrons located along the 
line of sight 

Xe = D^n. (7) 

The appropriate measure of integration is g{D) dD, which represents the probability that a CMB photon observed 
today scattered in the interval dD and has traveled freely since then. Here, g{D) is the visibility function and is 
related to the optical depth t{D) according to 

.(D)^-%e-"-K (8) 

The resulting polarization is [321 [3^] 

/6 /■ ^ 

(Q±zC/)(n) = ^ / dDg{D) ^ T^^iDn) ±2^2™ (n) , (9) 

m=-2 

where Q, U are the Stokes parameters of the CMB polarization and ±2^2171(11) are the spin- weighted spherical har- 
monics. These may be alternatively expressed in terms of the E and B polarization moments 

ztiBim) ±2Yim{n)- (10) 



IV. BUBBLE COLLISION POLARIZATION 



Armed with (|9]), we need rdc,e//(xe + i?efie) in order to find the polarization resulting from a bubble collision. This 
follows immediately from ([5]): 

r.c,e//(Xe+ i?efle) ^ ^ ^ ^^^^^^ _ ^^^^^^^^^ _ ^^^^ 
(L + ZdcjJ-O 
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FIG. 2. The geometry of the Earth's past hght cone projected on the xy plane of the decouphng surface, showing the primary 
reionization and decoupling scattering surfaces at 2 ~ 10 and z ~ 1050 respectively, and the decoupling sky as seen by an 
electron at some angle on those surfaces. The w direction is suppressed. 



Here — cos^e = — Xe)/Re is the cosine of the angle on the electron's CMB sky, and iie.c — {xc — Xe)/Re is the 
cosine of the angular radius of the collision disk on the electron's CMB sky (c./. Q)!^ 



The effective temperature anisotropy (11) is ready for substitution into the expression for the Stokes parameters. 
The planar symmetry of the temperature perturbation Tdc means that the only non-vanishing component of the 
quadrupole moment is 

2 




(12) 

where in the first equality we have used i?eMe,c = Xc — Xe — D^cfJ-c — D^fJ-, and a = Dg/Ddc is the ratio of the distance 
to the electron to the distance to the decoupling surface. The Stokes parameter U vanishes in these coordinates, again 
due to the planar symmetry, as do all the spin- weighted spherical moments of Q except those with m = 0. Therefore 
Bim = and the polarization is purely E'-mode — as expected for a scalar perturbation. Using 



^2o(Ai) = T\/^(l-A*'), (13) 



we find 



±2 



Qi(^) = i^{^] i'daD,,g{aDac){l-a)n-[^^^-^] }e{...), (14) 









""7 
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) r 

max 






^ If fJ,e,c > 1 the electron does not see the collision, if /Xe,c < —1 the electron's entire sky is taken up by it. 
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where the argument of the step function is the quantity in { }, and {AT/To)^^^^^ = X{Ddc — Xc) is the temperature 
perturbation at the center of the colhsion disk. For a hot spot [AT /Ti^)^^^ > and Q < (recaU that g = 
-{dr / dD)e^'' < 0), while for a cold spot Q > 0. 



V. ANALYSIS OF Q{e) 



With formula (14 1 in hand, our goal is to compute Q{9). In (14 1, Q{0) is determined up to the behavior oi g{aDdc), 
the visibility function in ([s]) that encodes the information about when CMB photons last scattered. It is normalized 
by 



/ g{D)dD = - f daDdcg{aDdc)^l. 



(15) 



Photons are most likely to scatter at roughly two times: within one scattering length of decoupling (at a redshift of 
approximately z ^ 1050), and at reionization {z ~ 10). 

The atomic physics at decoupling is well understood, and we can compute the visibility function there to good 
accuracy as a function of the cosmological parameters using [30] ■ However, the approximations we used in computing 
the temperature anisotropy from the collision break down at angular scales of about a degree and below, due to effects 
such as acoustic oscillations in the plasma between reheating and decoupling. Silk damping, the finite thickness of 
the decoupling surface, etc. We are currently engaged in studying these effects in detail |25) . For now, one must 
bear in mind that the signatures presented here receive 0(1) corrections on degree scales and below {£ > 200). This 
is particularly relevant for the contribution to polarization due to scattering at decoupling, since its details depend 
on the sub-degree structure of the temperature distribution. Nevertheless, the rough (> 1° scale) features of the 
signature we calculate are robust to these corrections. 

The physics responsible for reionization is less well understood, and our results are somewhat sensitive to the 
reionization model. Therefore we will proceed as far as possible with a model-independent analysis, and then study 
some representative reionization models. Using some reasonable assumptions we will see that we can make good 
estimates for Q{d) and quantify the expected effects from a collision into a unique signal that can be correlated with 
the effect on the CMB temperature map from the collision, e.g. the WMAP cold spot 



A. The Model-Independent Geometric Factor 



To isolate the dependence on the visibility function we rewrite ( 14 ) as 



Q{9)^ f daDdcg{Ddca) 
Jo 



(16) 



64 Uo^^ l-A^,^ '] 



1 — a 



©(■ 



= K{a,e) 

such that K{a,6) is dimensionless and contains all the recombination and reionization model- independent factors. 



1. Analysis of K( a, 9) 

The function K{a, 6) encodes the contribution of the collision perturbation to the Stokes paramater Q{9) due to 
Thomson scattering by a surface of electrons in our past lightcone at a comoving distance aDdc- This contribution 
is weighted by the visibility function to give the total Q, which is dependent on the model for recombination and 
reionization. It is thus instructive to examine the behavior of K in more detail. 

Thomson scattering at aDdc will lead to an annulus of Q-mode polarization on the sky, centered at = (the center 
of the spot). The inner radius of the annulus may be either zero (making it a disk) or non-zero, and its outer radius 
is always greater than 9^. By symmetry one can see that the polarization vanishes at = 0, which is encapsulated by 
the geometric prefactor (1 — /.t^). 
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The radii of the annulus are determined by the step function in (16). The argument of the step function is non- 
vanishing when 

Hcl-a iicl-a 

< fi< 1 . (17) 

a a a a 



For there to be a annulus rather than a disk, the upper bound in (17) must be <1. Physically, this corresponds to a 
region on the CMB sky for which all points on the electron's decoupling sphere are inside the collision region x > x^- 
In this case the electron sees a pure dipole on its sky, so there is no quadrupole moment to contribute to polarization. 
In Fig. [2] this happens when the circles corresponding to the electrons' past light cones are small enough, and the 
angle close enough to = 0, that they fit entirely in the region to the right oi x = Xc- In terms of the inequality in 



(17) this occurs when 

/Ltc < 2a - 1 . (18) 



When this is satisfied the inner radius is 



and in all cases the outer radius is 



cos"^ { — + ^ — " ) < Oc, (19) 



a a 



a 



cos-i f ^ - — ) > e,, (20) 



where the inequalities are saturated only if a = 1 (the scattering electron is at the decoupling surface) . 



B. Analysis of the Stokes Parameter 



We have analyzed the model-indepedent geometric factor involved in computing to the Stokes parameter. In this 
section, we will employ some fiducial models for the ionization fraction from reionization (combined with the well- 
understood ionization fraction from recombination) to compute the resulting Stokes parameter. The first step is to 
compute the visibility function. 



1. The Visibility Function 

Recall that the visibility function as a function of rcdshift is 

,(.) = --e-^(^) , (21) 

where t[z) is the optical depth, and g{z) is normalized so that — dz g{z) = 1. The optical depth to scattering out 
to cosmic time t is 

/to 
aTne{t)dt, (22) 

where ctt is the Thomson cross-section, rig is the number density of free electrons and <o is the age of the Universe. 
In terms of redshift this becomes 

r(z) = nl^ar ^^^^Xe{z')dz' , (23) 

where 71^°^ is the total number density of protons today (both boimd and free), Xe = ^i-e/n^pa is the ionization 

fraction, and H{z) = HQyjQ,r{l + zY + rim(l + z)'^ + ^k{^ + zY + fiA is the Hubble rate at redshift z. Here, Hq 
is the present Hubble rate, and the VLx (for x £ {r, m, fc,A}) are the fractional energy density of the Universe in 
radiation, matter, curvature and dark energy (cosmological constant) respectively. We assume a standard helium 
fraction of Yue — 0.24 and use the best fit values from [33] as reference cosmological parameters. Thus, our only 
other input is the ionization fraction, which characterizes the physics of decoupling and reionization. For a given 



9 




FIG. 3. Polarization patterns around cold (Q > 0) and hot {Q < 0) spots. 

ionization history it is straightforward to compute the visibility function. It is physically illuminating to trade z for 
the comoving distance D as the independent parameter in the visibility function using g{D)dD = g{z)dz where 

- f ^ = ^f dafia) , (24) 

and /(a) = ['nr+flma+^kO-'^+^AO.'^]~^^^ ■ The physics of decoupling is well understood, and we compute the ionization 
fraction using |30j . For reionization, we will take two representative models, the single reionization paramatrization 
used by CAMB [3S] , and a more complex double reionization model to qualitatively span the likely space of possible 
reionization histories 1361 . 



2. Single Reionization 



We take the ionization fraction at late times to have the form 



Xreiz) 



1 — tanh 



U — u(Zre) 



uiz) = {l + z) 



3/2 



(25) 



Vl + ZreA^ , 



where h, z^e and are input parameters. We set h = 1.08 in agreement with CAMB (taking the first reionization of 
helium into account) |35j . We choose values that produce a realistic ionization fraction at reionization [35 , normalize 
such that Tre = 0.09, Zre = 10.8 are in agreement with WMAP-7 year data, and choose = 0.5. The total ionization 
fraction is Xe{z) = Xre{z) + Xdc{z) where Xdc is taken from RECFAST [3D] with best fit values from the WMAP-7 
year data [34j . We plot the ionization history for this choice in Fig. |4j Equipped with this, we compute the visibility 
function g(aDdc) using the steps outlined above (taking D^c — D{z ~ 2000) as our reference distance scale), and then 
use ( [l4| to compute the Stokes parameter Q. 

We compute Q for a small spot with parameters similar to the WMAP cold spot (|Ar/To|max = 2 x 10~^ and 
an angular radius of 10°), and for a larger spot (| Ar/To|niax| = 5 x 10^^ and angular radius of 60°). We show the 
results for |Q(6')| in Fig. [5j The polarization patterns are radial for cold spots {Q > 0), and tangential for hot spots 
(Q < 0), see Fig. |3] For the small spot, the dominant contribution comes from scattering at decoupling and so is not 
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FIG. 4. The ionization fraction in the single reionization model with Zr^ — 10.8, Az — 0.5 and h — 1.08. We normalize such 
that Tre — 0.09; the early time ionization fraction is computed using RECFAST with best fit values from WMAP-7 year data 
|34| . Figure (a) is for all times and (b) zooms in on the reionization epoch. 
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FIG. 5. The Stokes parameter |(5(S)| vs. 6 (solid black line) for the single reionization model of Fig. |4] Figure (a) is the 
small spot I at/To I max — 2 x 10"'* and an angular radius of 10°. The dominant contribution for the small spot comes from 
scattering at decoupling. Figure (b) is the large spot with jAT/To|max = 5 x 10""" and an angular radius of 60°. For the large 
spot the sharp feature is the contribution from scattering at decoupling, the rest of the function comes from scattering after 
reionization. Q < for a hot spot and Q > for a cold spot, see Fig. |3] On each we display the estimated sensitivity for both 
the SPIDER (red dotted line) and Planck (blue dashed line) experiments, averaged over an annulus with inner angular radius 
8 and a width of 2°. 



sensitive to the physics at reionization. In contrast to this, for the large spot the dominant contribution comes from 
scattering at reionization, while scattering at decoupling leads to the narrow spikej^ Thus larger spots are sensitive to 
the physics at reionization. The Stokes parameters are robust under reasonable variations of the reionization history 
(within the single reionization model). 

To give the reader an estimate of the detectability of our signal we estimate the sensitivity level for two polarization 
experiments, SPIDER and Planck. We estimate the sensitivity by computing the number of pixels in an annulus with 
inner radius given by 6 in Fig. [sjja) and (b) and an angular width of 2°. On each figure for the Stokes parameter 
we display the resulting sensitivity level based on current estimates for SPIDER (red dotted line) in the 145 GHz 
band [57) and Planck (blue dashed line) in the 143 GHz band [3S]. For SPIDER we assumed a 25 day flight and 
and coverage of roughly 50% of the sky (which includes the area of the cold spot). For Planck the data is over the 
full expected 15 month lifetime of the instrument. For both the small and large spot the signal should be detectable 
based on current estimates of noise levels and detector sensitivity. 



We emphasize again that sub-degree features in our analysis are subject to 0(1) corrections, because we have not taken the relevant 
sub-horizon physics into account. We expect this spike to broaden to roughly a degree, and potentially to develop substructure |25) . 
However, the qualitative feature — a sharp spike at the edge of the temperature disk — should remain. 
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FIG. 6. Ionization fraction in the double reionization model with h/h = 1.1, zi = 4.6, Az = 0.5, Z2 = 20 and (72 = 2.9. We 
normalize such that h — 1.08 and Tre = 0.09 in agreement with WMAP-7 year data. The early time ionization fraction is 
computed using RECFAST with best fit values from WMAP-7 year data [53]. Figure (a) is for all times and (b) zooms in on 
the reionization epoch. 




FIG. 7. The Stokes parameter ](3(S)| vs. 9 for the double reionization model of Fig. |6] Figure (a) is the small spot 

1 at/To I max — 

2 X 10"* and an angular radius of 10°. The dominant contribution for the small spot comes from scattering 
at decoupling. Figure (b) is the large spot with |AT/ro|max = 5 x 10~^ and an angular radius of 60°. For the large spot the 
sharp feature is the contribution from scattering at decoupling, the rest of the function comes from scattering after reionization. 
Q < for a hot spot and Q > for a cold spot, see Fig. [3] On each we display the estimated sensitivity for both the SPIDER 
(red dotted line) and Planck (blue dashed line) experiments, averaged over an annulus with inner angular radius 9 and a width 
of 2°. 



3. Double Reionization 



We model double reionization using 



Xre{z) = 



1 — tanh 



u — u{zi) 



(26) 



where h, Zre, ^z, h, zi and ui are input parameters. We again normalize using h = 1.08. The other parameters 
are free, but must be chosen such that r^e = 0.09. As a typical example, we choose h/h = 1.1, zi = 4.6, A2 = 0.5, 
Z2 = 20 and 02 = 2.9. The ionization fraction is shown in Fig. [6j This function is a good approximation to that of 
|36j . We compute the Stokes parameter Q using the same spot parameters as before, and display this in Fig. [7j 

As in the single reionization case, we see that for small spots the dominant contribution is from scattering at 
decoupling, while for large spots it is scattering at reionization. The Stokes parameters are again robust under 
reasonable variations of the reionization parameters within the double reionization model. For the large spot, Q 
differs slightly from the single reionization model. We again display the expected sensitivities of SPIDER in the 143 
GHz band (red dotted line) and Planck in the 145 GHz band (blue dashed line) on each figure. 

The polarization pattern for the models plotted in Fig. [5] and Fig. [7] would be detectable by both Planck and 
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SPIDER given their instrumental sensitivities and angular resolutions. A smaller-scale polarization survey, with 
better resolution and sensitivity and focused on the cold spot, could improve on this. 

VI. DIFFERENTIABILITY FROM PRIMORDIAL FLUCTUATIONS 

The temperature and polarization anisotropies generated in the CMB from a bubble collision reflect an underlying 
planar symmetry arising from the direction picked out by the collision. This planar symmetry differentiates the 
polarization pattern generated via a bubble collision from those due to generic random Gaussian fluctuations. This 
is qualitatively illustrated in Fig. [8] where the temperature and polarization anisotropies from a bubble collision are 
contrasted with those from a Gaussian random field constrained to produce a temperature anisotropy identical to 
that from a bubble collision and a polarization pattern consistent with the angular power spectra Cf"^ and Cf^. 
For a Gaussian random field it is not necessary for the polarization pattern to be circularly-symmetric even if the 
temperature pattern is, and this reflects the fact that the projection from three-dimensional curvature perturbations 
to two-dimensional temperature perturbations is not unique — many three-dimensional patterns can produce the 
same two-dimensional pattern. Randomly generating this planar symmetry would be highly unusual, and so it is 
clear that a highly correlated and symmetric temperature and polarization pattern represents a non-trivial test of a 
bubble-collision origin for features in the CMB. We leave the important issue of quantifying the expected signal-to- 
noise of these correlated polarization features and a detailed analysis of detectability versus random fluctuations to 
future work [25] . 

VII. CONCLUSIONS 

In this paper we have computed the polarization of CMB photons due to a cosmic bubble collision. Together with 
the results of [5|, these results constitute a specific, quantitative prediction for the effects of a cosmic bubble collision 
on the CMB. However, there are other models that could produce signals of this type. Chief among them is of course 
concordance inflationary cosmology, which predicts the existence of random Gaussian fluctuations in the temperature 
map, along with their associated polarizations. Can one distinguish the effects of a bubble collision from those of an 
unusually large primordial random fluctuation? 

To answer this question in detail would require an analysis beyond the scope of this paper. However there is a 
sharp point to be made: the perturbation produced by the collision has a very characteristic geometry. It has planar 
symmetry, is zero outside the surface of the lightcone of the collision, and at least near the surface increases linearly 
with distance inside]^ The CMB temperature map provides a probe of a two dimensional slice of this three dimensional 
perturbation at decoupling. As we have seen in this paper, _E-mode polarization originates from and provides a probe 
of a different slice (or rather slices) . The information provided by polarization is independent — at least partially — of 
the information provided by temperature, and therefore the two together form a powerful check on the model. 

Imagine that a temperature disk of the type predicted by [S] were to be observed in the CMB (and in fact a disk 
with at least roughly the correct profile does indeed appear to be present [TM25] V Certainly such a feature could 
have been created by a random Gaussian perturbation — although perhaps with very low probability. But with what 
probability would such a random fluctuation produce both the correct temperature map and the correct polarization 
signal to correspond to the one computed here? The answer depends on many factors, among them uncertainties 
affecting the visibility function and reionization history of the Universe. But it is clear that at least in principle, one 
can have considerable leverage with which to rule out a random origin for such a signal. 

Another model that can produce cold or hot spots on the CMB temperature map are cosmic textures [35] ■ Because 
textures are late decaying, they do not produce a significant polarization signal |40j , and therefore can be distinguished 
from the collision model. Similarly if the temperature perturbation is an integrated Sachs- Wolfe effect due to a void 
(or overdensity) at 2; < 5 or so, its associated i?-mode polarization signal will be minimal. 

Of course the CMB is not the only tool available to study perturbations in the early Universe. Other observations 
sensitive to the primordial power spectrum, such as the Lyman-alpha forest [H], 21-cm radiation [121 US], and large- 
scale structure surveys (such as RH]) are also sensitive to the effects of bubble collisions [5l[6]. Because of the special 
characteristics of the collision geometry, these probes should produce results that are correlated with the CMB in a 
specific way. Since these measures access a three-dimensional volume, at least in principle far more data and far more 
statistical power remains to be tapped. If indeed there is a feature in the CMB sky due to a bubble collision that is 



We remind the reader again that this statement is vaUd only on angular scales larger than about a degree. On smaller scales a more 
detailed analysis is required 1251 . 
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FIG. 8. The upper panels are (from left to right) the circularly-symmetric T, Q, and U patterns resulting from the three- 
dimensional planar symmetry of a bubble-collision generated cold spot with an angular radius of 10 degrees and a peak 
temperature deficit AT/T — —2 x 10^*. The lower left panel is an identical circularly-symmetric T pattern but the Q and U 
panels are now a realization of a Gaussian random field statistically consistent with the T pattern, and the polarization power 
spectra Cp^ and Cf^ . This illustrates that generic polarization patterns consistent with a two-dimensional cold spot need not 
have the full three-dimensional planar symmetry expected from a bubble collision. 

sufficiently bright to be observed, we are optimistic that a combination of these data would suffice to prove its nature 
beyond reasonable doubt. 

We close with one final comment. In this paper we have focused on the effects of a single collision. But the results 
for the probability and distribution of characteristics of bubble collisions derived in [7J had an interesting feature — the 
measure on the size of the collision disks is flat in fic = cos 0c oc Xc- The reason for this is fairly obvious in retrospect: 
because of inflation, only a tiny region of the last scattering surface is visible to us today. Collision lightcones slice 
through this surface at various positions. While the distribution of their locations Xc,i is not uniform, because we can 
only see a very small part of the surface, it is uniform to a very good approximation over the range we can observe. 
Therefore disks with very small 9c (and by the same token, disks with very small tt — 9c) are unlikely. The cold spot 
in the CMB has 9c ~ 10° in angular radius. Since the measure is flat in cos6'c, the probability for this to be the only 
collision disk visible in our sky is small. Instead, it may simply be the brightest.... which means that a search for 
larger, fainter disks is well worth doing |45) . 
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